Spin squeezing via quantum feedback 
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We propose a quantum feedback scheme for producing deterministically reproducible spin squeez- 
ing. The results of a continuous nondemolition atom number measurement are fed back to control 
the quantum state of the sample. For large samples and strong cavity coupling, the squeezing param- 
eter minimum scales inversely with atom number, approaching the Heisenberg limit. Furthermore, 
ceasing the measurement and feedback when this minimum has been reached will leave the sample 
in the maximally squeezed spin state. 

PACS numbers: 42.50.Dv, 32.80.-t, 42.50.Lc, 42.50.Ct 



Squeezed spin systems Q of atoms and ions have at- 
tracted considerable attention in recent years due to the 
potential for practical applications, such as in the fields of 
quantum information [p| and high-precision spectroscopy 
||. Quantum correlations of squeezed spin states out- 
perform classical states in analogy with squeezed optical 
fields. Moreover, squeezed spin sates are also multiparti- 
cle entangled states [||. Recent proposals for their gen- 
eration include the absorption of squeezed light ||], col- 
lisional interactions in Bose- Einstein condensates j|, |), 
and direct coupling to an entangled state through inter- 
mediate states such as collective motional modes for ions 
or molecular states for atoms gj. 

Other proposals create spin squeezing via quantum 
nondemolition (QND) measurements B, [II], O. A strik- 
ing recent achievement of QND measurements is the 
entanglement of two macroscopic atomic samples [ p"2[ . 
These QND schemes produce conditional squeezed states 
that are dependent on the measurement record. On the 
other hand, unconditional squeezing would ensure that 
the state is deterministically reproducible. Maimer Jfj| 
has shown that alternating QND measurements and inco- 
herent feedback can produce sub-Poissonian number cor- 
relations. However, that work does not treat the quan- 
tum effects of the measurement back action or the feed- 
back on the mean spin (which is assumed to be zero). 
Hence it cannot predict the strength of the entanglement. 

In this Rapid Communication, we suggest achieving 
spin squeezing via feedback that is coherent and contin- 
uous. We consider a continuous QND measurement of 
the total population difference of an atomic sample. The 
results of the measurement, which conditionally squeeze 
the atomic sample, are used to drive the system into the 
desired, deterministic, squeezed spin state. This involves 
amplitude modulation of a radio-frequency (rf ) magnetic 
field, where the feedback strength varies in time such that 
the mean number difference is kept at zero. 

An ensemble of N two-level atoms can be described by 
a spin- J system RM, i.e., a collection of 2 J = N spin- 



i particles. The collective spin operators are given by 



J a — J2k=i ho~£'(a — x,y,z), where are the Pauli 
operators for each particle. Thus, J z represents half the 
total population difference. Coherent spin states (CSS) 
have variances normal to the mean spin direction equal to 
the standard quantum limit (SQL) of J/2. Introducing 
quantum correlations among the atoms reduces the vari- 
ance below the SQL in one direction at the expense of the 
other jlj . Such squeezed spin states can be characterized 
by the squeezing parameter B 



4 =AT(AJ ni ) 2 /(<^ 2 } 2 + (J n3 ) 2 ), 



(1) 



where rij(i = 1, 2, 3) are orthogonal unit vectors. Systems 
with £ 2 < 1 are spin squeezed in the direction n and also 
have multiparticle entanglement Q. 

Let the internal states, |1) and |2), of each atom be the 
degenerate magnetic sublevels of a J = | state, e.g., an 
alkali ground state. Each atom is prepared in an equal 
superposition of the two internal states, thus giving a 
CSS of length J in the cc-direction. The atomic sample is 
placed in a strongly driven, heavily damped, optical cav- 
ity, as shown in Fig. |](a) . The cavity field is assumed to 
be far off resonance with respect to transitions probing 
state 1 2), see Fig. ||(&). This dispersive interaction causes 
a phase shift of the cavity field proportional to the num- 
ber of atoms in |2). Thus, the QND measurement of 
J z (since N is conserved) is effected by the homodyne 
detection of the light exiting the cavity jl5| . 

This interaction is defined by the Hamiltonian hxJzb^b 
where b,b^ are the cavity field operators and x = <7 2 /4A, 
with one-photon Rabi frequency g, and optical detun- 
ing A For strong coherent driving we can use the 
semiclassical approximation b — > i\[3\ +6, where now h 
represents small quantum fluctuations around the classi- 
cal amplitude i\j3\. The interaction is thus 

H int = h x \(3\J z (-tb + tb^), (2) 

where we have chosen an initial J z splitting of — xl/^| 2 - 
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FIG. 1: (a) Schematic experimental configuration. A cav- 
ity field of amplitude i\j3\ interacts with the atomic sample. 
The current Ic(t) from the homodyne detection of the cavity 
output, damped at rate k, is combined with A(t) produced 
by a signal generator (SG). The combined signal controls the 
amplitude (AM) of an rf magnetic field that, together with a 
static field, drives J y . (6) Single atom diagram. The static 
B field lifts the degeneracy of the magnetic sublevels. The 
far-detuned cavity field i\/3\ monitors the collective popula- 
tion in state \2) (and hence J z ). The rf driving field, applied 
perpendicularly to the static field direction, induces magnetic 
dipole transitions between |1) and |2) (thus driving J y ). 




FIG. 2: Schematic (a) y — z and (b) x — y projections of 
the quasiprobability distributions for the spin state. The spin 
states are represented by ellipses on a sphere of radius J. 
The initial CSS, spin polarized in the x direction, is given 
by state 1. State 2 is one particular conditioned spin state 
after a measurement of J z , while state 3 is the corresponding 
ensemble average state. The effect of the feedback is shown 
by state 4. A rotation about the y axis shifts the conditioned 
state 2 back to (J z ) c = 0. The ensemble average of these 
conditioned states will then be similar to state 4. 



Following the procedure of Sec. VII in Ref. |L6[ we can 
adiabatically eliminate the cavity dynamics if the cavity 
decay rate k ^> x|/3|J 2 , which requires k 3> xl/Jjv'iV 
(since the initial AJ Z = \J J/2). The evolution of the 
atomic system due the measurement is thus 



i> = MV[J z ]p, 



(3) 



where M = 4x 2 \[3\ 2 /k is the measurement strength 
[equivalent to 2D in Eq. (22) of Ref. (XB)], and V\r\p = 
rpr^ — (r^rp + pr^r)/2. This equation represents deco- 
herence of the atomic system due to photon number fluc- 
tuations in the cavity field, with the result of increased 
noise in the spin components normal to J z . 

The effect of Eq. @ on the cavity field is a phase 
shift proportional to J z , and thus the output Homodyne 
photocurrent is given by |16) 



I c (t) = 2VM(J Z 



(4) 



where ((t) is a white-noise term satisfying E[£(i)£(i')] = 
5(t — t') and E is the ensemble average. The conditional 
master equation for the atomic system is then Jlq] 



(5) 



dp c = dtMV[J z ]p c + VMdW(t)H[J z ]p c , 



where dW(t) — £(t)dt is an infinitesimal Wiener incre- 
ment and H.[r]p = rp + pr^ — Tr[(r + r^)p]p. 

The effects of this evolution on the initial CSS are a de- 
crease in the variance of J z with corresponding increases 
for J y and J x (i.e., spin squeezing), as well as a stochastic 
shift of the mean J z away from its initial value of zero. 
This shift, indicated by state 2 of Fig. ||, is equal to 



d(J z ) c = 2VMdW(t)(AJ z ) 2 c m 2VMI c (t)dt(J 2 ) c . (6) 



Here the approximation assumes that (J z ) c = 0, which 
will be relevant when feedback is in place. 

The average or unconditioned evolution, Eq. (|^), is 
simply recovered by taking the ensemble average of all 
possible conditioned states, i.e., p(t) = ~E\p c (t)]. This 
leads to a spin state with J z variance equal to J/2. In 
other words, the unmonitored measurement does not af- 
fect J z and the squeezed character of the individual con- 
ditioned states is lost, indicated by state 3 in Fig. |2[ 

To retain the reduced fluctuations of J z in the average 
evolution, wc employ a feedback mechanism that uses 
the measurement record to continuously drive the sys- 
tem into the (J z ) c — squeezed state. The idea is to 
cancel the stochastic shift of (J z ) c due to the measure- 
ment. This simply requires a rotation of the mean spin 
about the y axis equal and opposite to that caused by 
Eq. (H), as illustrated in Fig. |^. The feedback Hamilto- 
nian must therefore take the form 



H {b (t) = X(t)I c (t)J y /VM, 



(7) 



where X(t) is a time- varying feedback strength. This feed- 
back driving can be implemented by modulating an ap- 
plied rf magnetic field [[l7|, as shown in Fig. |l|. 

Following again the methods of Ref. |l6| to find the 
total stochastic master equation, we can calculate the 
conditioned shift of the mean J z due to the feedback. 
Again using the assumption that (J z ) c = we have 



d(J z 



fb 



-X(t)I c (t)dt(J x ) c /VM. 



(8) 



Since the idea is to produce (J z ) c = via the feedback, 
the approximations above and in Eq. (^) apply and we 
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find that the required feedback strength is 
\(t) = 2M(j 2 z ) c /(J x ) c . 



(9) 



This type of feedback control is essentially a form of state- 
estimation-based feedback |lgj |. Although Eq. (|?]) looks 
like direct current feedback, the strength of this feed- 
back (^) is determined by conditioned state expectation 
values. I c (t) only appears directly in H{b due to the as- 
sumption that the feedback works and so ( J z ) c = 0. 

Being dependent on conditioned expectation values, 
which are computationally very expensive, Eq. (^) is not 
practical in an experimental sense. What is required is a 
predetermined series of data points or ideally an equation 
for X(t), like in Fig. [p. To find a suitable expression we 
begin by assuming the feedback is successful and replace 
the conditioned averages by ensemble averages, 



\(t)~2M(j*)/(J x ) 



(10) 



This approximation will be valid if the unconditioned 
state has high purity since then it must comprise of nearly 
identical highly pure conditioned states. 

The evaluation of both the purity (Tr[p 2 ]) and the av- 
erages in Eq. ( |io|) ((A) = Tr\Ap\) requires the uncondi- 
tioned master equation (ME) [[l6[ 

p = MV[J z ]p-i\(t)[J y ,J z p + pJ z ] + ^V[J v ]p. (11) 

The terms in this equation describe, respectively, the 
noise due to the measurement back-action, the feedback 
optical driving, and the noise introduced by the feed- 
back. The state determined by Eq. (11), with X(t) given 
by Eq. (|Io|), has a purity very close to one (see below). 
Since the state is very close to a pure state, we are justi- 
fied in applying the approximation of Eq. (|To|). 

Note that Eq. ( |TT| ) describes the exact unconditioned 
evolution of the atomic system where the feedback 
strength is arbitrarily defined by A(i). Equation jlO| ) 
thus describes one particular feedback scheme, however, 
it can only (easily) be evaluated numerically. To find 
an approximate analytical expression we look at when 
(J) ~ O(J), for which the atomic sample remains near 
the minimum uncertainty state ( J 2 ) ( J 2 ) « J 2 /4. This 
is equivalent to a linear approximation represented by 
replacing J x with J in the commutator [J y ,J z ] — iJ x , 
which allows us to calculate (Jy) directly from the ME 
(and hence (J 2 ))- The decrease of (J x ) from J is then 
related to the increase of {Jy) from J/2 [see Fig. |^(6)] 
due to the measurement back action. Using these ap- 
proximations we obtain 

A(i) re Me Mt/2 (l + 2JMt)-\ (12) 

We can analytically approximate the degree of squeez- 
ing produced by the particular feedback scheme repre- 
sented by Eq. (|l2|). For our model Eq. ([!]) becomes 

N{jf) \{t)J 



e Mt {l + 2JMty 



M(,L 



(13) 



This leads to a minimum at t* pa 1 /M of 



(14) 



Thus, the minimum attainable squeezing parameter 
asymptotically approaches an inverse dependence on the 
sample size, i.e., the Heisenberg limit figf . 

The approximations leading up to Eqs. ( |l0| ) and (]l2]) 



can be justified by numerically solving the ME (11) for 
a given A(t), for example, using the Matlab quantum 
optics toolbox |po|| . The approximation of Eq. (10) can 
be tested by calculating the purity for p described by 
the ME with this particular feedback. The expectation 
values in A(t) are found by iteratively solving the ME 
[updating X(t) each time step], and thus we also have 



e z = 2j{Jl)/{j x y 



(15) 



The results of this simulation are shown in Fig. |^(a), 
where the purity is given by the dotted curve and £ 2 
is curve A. Clearly, the purity remains near unity for 
times of interest. This implies that the measurement and 
feedback scheme has worked to produce nearly identical, 
nearly pure, conditioned states [for times 0(M~ 1 )}, and 
we are therefore justified in using Eq. (|Io|). 

The further approximations to obtain the analytical 
Eq. ( [l2| ) are also good, as shown by curve B in Fig. ||(a), 
where the fit to curve A for the early evolution and min- 
imum is nearly perfect. We are not interested in later 
times since the idea is to cease the measurement and 
feedback when the minimum £ 2 is reached. The analyt- 
ical expression for the squeezing parameter, Eq. (plf), is 
also plotted as curve C in Fig. ||(a). Although the min- 
imum is not a perfect fit to the exact numerical results, 
it has the correct order of magnitude and so we expect 
the J -1 scaling of Eq. ( |l~4| ) to be correct. 

The scaling of £^ in is obtained numerically from solu- 
tions of the ME [and thus Eq. (|l|)] with feedback de- 
scribed by Eq. Jl^). This A(t), shown above to be a good 
approximation, is the suitable form for experimental real- 
ization. The numerical results, along with the analytical 
expression ( [Ti] ) , are plotted in Fig. ||(&) which clearly veri- 
fies the J -1 dependence. The analytical coefficient (e/2) 
represents an error of ~ 18% compared to the numeri- 
cal fit, and the optimum time (i») is also slightly out as 
shown by Fig. [|(a). Nevertheless, these errors only apply 
to scaling coefficients, not to the scalings themselves. 

Experimentally, the limit to squeezing will be domi- 
nated by spontaneous losses due to absorption of QND 
probe light. The rate of this loss is Njtl 2 /4A 2 , where 
7 is the spontaneous emission rate and fi = g\(3\. To 
reach the Heisenberg limit (requiring a time = 1/M) 
we want the total loss to be negligible, i.e., we require 



A 4 X 2 



1 



g 2 - N a?7, (16) 



which is the very strong coupling regime of cavity QED. 
Note we also required n ^> x|/3| V^V for the adiabatic 
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FIG. 3: (a) Time dependence of the purity = Tr[p 2 ] (dashed curve), and the squeezing parameter £| (curves A, B, and C). The 
purity and curve A are the results for a feedback scheme defined by Eq. (|l^) . For comparison, Eq. (^5|) is also simulated for the 
scheme defined by Eq. (h2|) and the result is curve B. Finally, curve C is the approximate analytical expression for squeezing 
parameter given by Eq. (O). In all cases J = 25. (b) J dependence of the squeezing parameter minimum, £min- Plotted are the 
results of numerical solutions of the ME, and hence Eq. (fig), with X(t) given by Eq. ( |l2| ) (circles) and X(t) x 120% (diamonds). 
These values approach 1.665/J and 1.744/J, respectively. The dotted line is the lower bound for single-shot feedback also with 
a 20% error, while the analytical result for perfect continuous feedback, Eq. dlj), is the solid line. 



elimination of the cavity field, and to satisfy both we thus 
require A ^> 7 TV 3 / 2 . It is not surprising that Eq. ( |l6| ) is 
the same requirement as for Andre and Lukin's model 
implementing the countertwisting Hamiltonian [Q, 
since writing Eq. ( |TT| ) in Lindblad form reveals such a 
term. Similarly, the condition for achieving some squeez- 
ing, i.e., £ 2 < 1, will be g 2 N > ivy. Further, we have cal- 
culated that a free space model [also given by Eq. ( pd| ) 
but with M equal to N in Ref. will also produce 
some squeezing, although the Heisenberg limit cannot be 
reached since by i» = 1/M all atoms will be lost from 
the sample. 

Figure also indicates that our continuous scheme 
is very robust to any experimental errors in the feedback 
strength, as opposed to a single-shot method. The latter 
approach consists of a single (integrated) measurement 
pulse (see e.g., Ref. followed by a single feedback 

pulse. If there is a relative error of e in the feedback 



strength, this will induce an error term (A J| rr ) 



''J/2, 



which will dominate the total variance for J ^> 1. Thus 
£min wu l have a lower bound of e 2 , and will never be bet- 
ter than J -1 . On the other hand, as shown in Fig. ||(6), a 
large (20%) error in X(t) for continuous feedback does not 
affect the J -1 scaling. We have also found this theoreti- 
cal scaling to be unaffected by inefficient measurements. 
Finite feedback delay time will also have a limited effect 
as long as it is faster than (JM) _1 . 

This Rapid Communication has presented a scheme 
for producing a spin squeezed atomic sample via QND 
measurement and feedback. The advantage over previ- 
ous QND schemes ||, 10 is that it provides uncondi- 
tional, or deterministically reproducible, squeezing. For 
very strong cavity coupling, the theoretical squeezing ap- 
proaches the Heisenberg limit £ 2 ~ 1/JV, while some 
squeezing will be produced at weaker coupling and even 
in free space (thus presenting a simple experimental test 
for quantum feedback). This indicates a stronger squeez- 



ing mechanism than collisional interactions in a Bosc- 
Einstein condensate where the scaling is N~ 2 / 3 j|, [|. 
Furthermore, by ceasing the measurement when this min- 
imum is reached, the maximally squeezed state could be 
maintained indefinitely. 

L.K.T and H.M.W wish to acknowledge inspiring dis- 
cussions with Klaus M0lmer and Eugene Polzik. 
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